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Abstract 


The  IR  laser  excitation  of  an  adbond  la  studied  by  both  a  quantum  mechanical  generalised 
■aster  equation  approach  and  a  classical  generalised  Langevia  approach.  The  role  of  phonons 
in  the  energy  flow  between  the  adbond  and  the  surface  is  considered.  The  latter  approach 
looks  further  at  local  heating  via  direct  excitation  of  surface  atoms.  It  is  seen  that  the 
Markovian  approximation  is  in  general  inadequate,  and  that  local  heating  is  an  important 
mechanism  for  desorption. 


Introduction 


The  coherence,  monochromaticity  and  high  energy  density  of  lasers  have  recently  been  ex¬ 
ploited  to  Influence  dynamical  processes  occurring  at  a  solid  surface.  The  many  diverse 
applications  include  catalysis,  localised  melting,  chemical  vapor  deposition  and  charged 
particle  generation.  While  same  of  the  macroscopic  aspects  of  the  phenomena  occurring 
during  the  laser-stimulated  surface  processes  have  been  studied  theoretically,  a  fundamental 
microscopic  understanding  is  lacking.  In  this  paper  we  shall  address  the  role  of  phonons  in 
the  excitation  and  relaxation  of  adspeeles.  In  Section  II  we  present  a  quantum  mechanical 
approach  utilising  the  generalised  master  equation  to  obtain  the  probability  destribution  of 
a  laser-excited  adbond.  In  Section  III  we  present  a  classical  generalised  Langevin  approach, 
with  application  to  flash  desorption.  Bern,  laser  excitation  of  surfaoe  atoms  in  addition 
to  the  adbond  is  also  considered.  Conclusions  are  presented  in  Section  IV. 


Qusntwi  Generalised  Master  Equation 

In  a  recent  paperl  we  formulated  s  general  theory  of  vibrational  excitation  of  an  adatom- 
surface  bond  (adbond)  by  an  IB  laser,  accompanied  by  relaxation  via  vibrational  energy  ex¬ 
change  between  the  adbond  and  lattice  vibrations.  That  formalism  is  based  on  a  first- 
principles  Hamiltonian  approach  utilising  first-order  perturbation  theory,  and  the  effective 
potential  appearing  in  the  sero-order  Hamiltonian  for  the  adbond  is  assumed  to  be  a  sum  of 
pair  potentials  between  the  adatom  and  the  lattice  atoms  at  their  equilibrium  positions. 

For  the  sake  of  simplicity,  we  modeled  a  process  in  which  a  short  pulse  of  radiation  excites 
the  vibrational  states  of  the  adbond  which  subsequently  decays  via  energy  transfer  to  the 
phonons.  Bore  we  extend  that  earlier  work  to  include  more  general  types  of  laser  excitation 
and  introduce  an  alternative  type  of  effective  potential. 

Time  Evolution  of  the  -A^ond.a«n<i 

Be  are  primarily  Interested  in  the  time  evolution  of  the  vibrational  states  of  the  adbond, 
which  la  excited  by  IR  laser  radiation  and  is  coupled  to  the  phonon  modes  of  the  solid.  Ve 
choose  a  laser  whose  frequency  is  close  to  a  normal  mode  frequency  of  the  adbond  but  sot  to 
any  of  those  of  the  solid.  The  total  Eamlltoalaa  of  the  system  can  be  written  as 

Mere  T-  is  the  kinetic  energy  operator  for  the  adpertlcle,  is  the  interaction  between 
the  adpertlcle  sad  the  A-th  lattice  stem,  T.  is  the  kinetic  Kergy  of  the  letter,  V,t*  is 
the  interaction  energy  between  the  A-th  and  A'-th  lattice  atoms,  I.  is  the  Hamiltonian  for 
the  free  radiation  field  aad  I.,  is  the  Interaction  between  the  adbond  and  the  laser  radl- 
atlen.  The  study  of  the  time  Point ioa  of  the  vibrational  states  of  tho  adbond  then  in¬ 
volves  extracting  the  relevant  information  from  tho  equations  of  motion  of  operators  0  tor 
the  whole  system 

g  -  -i[J,0]  -  -1L0  ,  (» 

Xa  addition,  wo  also  need  to  specify  a  sot  of  Initial  conditions  compatible  with  the  history 
of  tho  system  prior  to  tho  time  tho  experiment  is  initiated. 


/ 
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;  1  z'Ji 


We  rewrite  oar  Hamiltonian  in  the  form 


where 


- T.  *  F*‘  • 


■>  ■  i<T*  *  ’ 


•»  ■  [<».«-%«>  •  <«> 

T-t  la  the  thermodynamic  average  of  over  the  Initial  equilibrium  configuration  of  the 
pfionon  fieldt 

*at  *  <<Val>>p  *  TVp0”^  *  (7) 

Pn(0)  being  the  Initial  equilibrium  density  operator  for  the  phonon  field.  The  adbond  is 
now  described  by  the  Hamiltonian  Ha.  The  radiation  field  £s  assumed  to  be  that  due  to  a 
single-mode  laser  of  frequency  aq,  and  polarisation  rector  c.  Accordingly, 

Br  ■  du^ata  ,  (8) 

and,  in  the  dipole  approximation,  e±l**?  ■  1,  with  effeetlre  adatom  charge  eg, 

**s(a-a+>  ,  (»> 

where  a+,a  are  creation-annihilation  operators  for  photons,  t  is  a  quantisation  volume  and 
k  is  the  ware  rector  of  the  radiation  field. 

Using  the  projection  operator  technique,2  we  arrlre  at  a  set  of  coupled  lntegrodlffereatlal 
equations  for  the  matrix  elements  of  a  projected  density  operator  p  for  the  system: 


PSa(t) 


-  £  jdf  He  ^«A8g,<f)A88,»e1“88’tj>sa(t-t 

-  He|<<AM,(f)A8S.»e"1“S8,t  Jp8'8'<t-f) 


where 


Pu<t)  -  <8|p(t)|8> 


p<t)  -  TrpTrrW(t)  (18) 

with  TrB  and  Tr*  representing  traces  orsr  the  phonon  and  radiation  spaces.  |g>  is  an  eigen¬ 
state  of  with  energy  Ig ,  tugS  •  is  giren  as 

"gg.  -  <VV>/fi  »  (18> 

and  V(t)  is  the  density  operator  of  the  whole  system. 


in  most  oases  of  chemical  interest,  some  charge  transfer  takes  place  between  an  atom  and 
a  surface  during  adsorption.  Is  shall  assume  the  atom  to  bare  an  effective  charge  en  and 
the  charge  transferred  to  the  solid  to  be  distributed  throughout  the  latter  so  that  it  is 
still  essentially  neutral.  Image-charge  effects  will  be  important  for  the  case  of  metals, 
but  we  shall  think  in  terms  of  aoametals  here  to  keep  the  treatment  simple. 

The  wibratioaal  levels  of  the  adatom- surface  bead  are  dictated  by  the  effective  potential 
seen  by  the  adatom.  This  potential  is  quite  complex,  being  a  result  of  electrostatic  and 
exchange  forces  due  to  the  nuclei  and  electrons  of  the  solid,  and  is  not  known  for  any  real 
system.  A  simple  representation  of  tho  effective  potential,  however,  can  be  obtained  by 
needing  independent  pairwise  interactions  between  tho  adatom  and  each  lattice  atom  and 
summing  over  all  lattice  atoms.  The  advantage  of  such  a  procedure  is  that  tbs  influence  of 
lattice  atoms  more  than  a  few  lattice  constants  sway  from  the  adatom  in  negligible  (tee  lat¬ 
tice  atoms  should  be  sufficient  for  most  solids).  Xf  the  sms  of  pair  potentials  oaa  farther 


:  1201 


b*  fitted  to  a  single  analytic  potential,  the  vibrational  levels  of  the  adatom- surface  bond 
reduce  to  those  of  a  diatonic  molecule. 

Ve  define  a  one-dimensional  coordinate  system  with  origin  at  the  mean  position  of  the 
•  outermost  lattice  atom.  Let  z  be  the  position  of  the  adatom  and  z.  the  position  of  the  i-th 
lattice  atom  measured  in  unite  of  [0(O)]-1  to  be  defined  later, 


*1  "  *t  *  m  ♦  «£  . 


1-1, a.. . 


«1  i»  the  displacement  of  the  1-th  lattice  atom  from  its  equilibrium  position  s(0)  • 
all  in  units  of  lfl(0)]-l#  assuming  a  monatomic  lattice  with  uniform  spacing  g.  The 
interaction  potential  at  the  position  of  the  adatoa  is  then  a  sub  of  pairwise  Morse  poten- 
tlals  between  it  and  all  the  atoms  of  the  chain, 

V<s,{*t>)  -  £  v(z-zt)  *  £  Vat  ,  (15) 

,  \  .<o,r-^0>^<— 1>  -{*-*J0>+(l-l)a-ut}1 

»(*-*!>•  D;;L*  -2e  0  I,  (16) 

where  z&0)is  the  position  of  the  minimum  of  »(z),  0<O)is  now  identified  as  the  Horse  ex¬ 
ponent  parameter  and  D£0>  Is  the  well  depth  of  the  Morse  potential.  The  zero-order  adbond 
Hamiltonian  becomes 

_  (g(Q)jt)3  r  -a(z-z0)  -(z-zo)] 

0  j.  <«) 

wh#r?n?e  sad  Sq  are  effective  parameters  which  involve  the  pair-potential  parameters  d£0) 
end *|r 'and  the  lattice  parameters  a  and  «ua»,  where  the  last  quantity  is  the  mean  square 
dlsplac  meant  of  a  lattice  atom. 


Within  this  model,  the  generalized  master  equation,  Bq.  (10),  becomes 


PM<*> 


■  j.  ?  {*•  J*ss: <■*■  •  > (|tTSs- 


where  the  A's  are  time- independent  coefficients t 


.»»•<„.>  - 1 [.-»<—40>>]  f -•<*- 4°» 

Alf'a*a  1  jiy1  j®*' 

x  exp  £-y  (  i-1  )a-u  •  (  i  •  -1  )a  ♦  H—tjj' .  «ua»j 


and  the  T's  are  phonon  dlepl 


it  correlation  fuzetloas, 


(t )  -  uu’«u1(t)uj,»  .  (ac 

i  term  p_  M(t)  is  the  contribution  to  &z(t)  due  to  the  radiation  field,  leplaoiag  b. 
P8,  we  SiPwrite  Bq.  (It)  as  ^ 

*sCt)  -JdfJ^I^  ^M,(t,)Pi,Ct-t,)-B«ia(t,)Pg(t-f^  ,  (M 


where  K  oonaists  of  contributions  dee  to  the  phonons  and  the  radiation  field, 

Xgg.Ct’)  -  t£2?(f)  ♦  lj5?(t'>  .  (SI 

For  n  system  such  as  oxygen  on  germanium,  the. amplitude  of  E&82  is  seen  in  Tig.  1  to  be 
ijr  orders  of  mngzltnde  larger  than  that  of  Kgfr  dee  to  a  1  l/m*  laser.  However,  the 


time  t  <pe> 


Fig.  1.  Kernel  functions  and  K&f *  for  0/0*.  On*-,  two-  sad  thrve-pbonon  procsss «s 

srs  Included.  The  laser-adodnd  ■!— stoh  Is  3  as*1. 


"frequency"  of  is  dictstsd  by  ths  Dsbys  frequency  of  tho  solid,  whereas  that  of  k£§? 

depends  on  the  decree  of  resonance  between  the  laser  aad  the  levels  of  the  adbond  and  oaa  be 
■ads  arbitrarily  swill  within  the  rotstlag-wave  approximation.  For  probability  functions  P8 
that  chance  slowly  durlnc  the  vibrational  period  of  typical  phonons,  the  contribution  due 
to  K»)  will  be  subject  to  oaneellatloae,  whereas  that  due  to  Kbit  will  not.  Whether  thla 
dynamical  effect  can  brides  the  cap  due  to  the  difference  la  nagaltude  of  the  amplitudes. of 
K$g/  aad  Kggf  depends  on  the  sharpness  of  the  laeer-adboad  resonance,  the  details  of  E&gf 
aad  the  total  Irradiation  tine.  -  • 

An  Important  consequence  of  the  nature  of  il  that  the  dynamics  Is  seen  to  be 

non-Mar kovlaa  since  there  Is  no  Indication  that  XA§/  In  localised  In  tine.  Is  ceneral,  one 
expects  very  complicated  sou local  effects  due  to  the  combined  laser  phonon  fields,  and  great 
caution  must  be  used  when  applying  a  Markovian  approximation  to  seek  processes.  Mumerlcal 
Investigation  of  the  detailed  time  evolution  of  a  number  of  laeer-adbond-eolid  systems  Is  Is 
progress  aad  should  provide  seme  mere  definitive  answers. 


In  thla  section,  we  shall  formulate  the  problem  of  flash  desorption  via  a  Langwvln  ap¬ 
proach.  We  model  the  problem  as  an  Incoherent  pro ones.  Thus,  the  motions  of  the  phonons 
nerve  as  a  beat  bath  for  the  adspeelen. 

WS  start  with  the  equation  of  motion  for  the  1-th  atom  la  a  three  dimensional  solid 


Vi  -  -J  Vi  ♦  '1 ♦  j  “u'5  ♦  '1 X  -1J*5  ♦  t5  X  ‘»J'J  *  I  •«*’  ♦•••  <*»> 


where  m.  aad  r.  are  the  nass  aad  position  of  the  1-th  solid  atom,  aad  kn,  aij,  bjj ,  etc. 
are  tharelastli  and  higher-order  nonlinear  fores  const ants,  for  simplicity,  let  us  consider 
the  atoms  to  have  unit  masses  aad  also  assmee  tbs  ssnlssal,  nsnllaccr  effects  to  he  un im¬ 
port  aet.  Moose,  It.  (MM)  reduces  to 


•»  *  ‘u*J  ♦  *1  }  ‘ll'j  *  'l  }  *ljr5  *  *5  J,  *U’ )  *— 


WS  rewrite  Mq.  (M)  in  matrix  nstatlsn  as 


(14) 


?■-!£♦  s**r  *  s'fty  ♦  x’gc  ♦••• 


Izdl 

(25) 


with  j'  and  j*  diagonal  matrices,  r£t  -  rt  and  v^  -  vt  -  rf . 

Vw  now  define  projection  operators  P  and-Q  each  that 

X<t)  -  P  r(t),  s(t)  -  (l-P)r(t)  -  Q  r(t),  PQ  -  0,  PP  -  P  (25) 

{y(t)>  are  coordinates  of  those  solid  atoas  Interacting  with  the  adspecles,  and  (s(t)>  are 
coordinates  of  the  remaining  solid  atoas.  Applying  these  projection  operators  to  Sq.  (25), 
n  obtain 

f(t)  -  Hfcp  ♦  j*(t)  jpp  ♦  j'(t)|pp  x'(t)  ♦  **-l3f(t)  ♦ 

I-fcpq  ♦  j*(t)  {pq  ♦  x'(t)  bpq  n*(t)  ♦  •••];(t>  (27) 

and 

S(t)  ■  (-kqq  ♦  g'(t)  sqq  *  §'(t)goQ  j'(t)  +  • « * ]s(t)  ♦ 

l-Jqp  +  f’(t)  sqp  +  f(t)^Qp  x*(t)  ♦  •••]x(t)  (28) 

*'(t)  and  y'(t)  are  diagonal  aatrlces  with  s{1(t)  -  a^t)  and  yit(t)  -  y,(t),  and  bp  -  P», 
|pq  “  PfcQ,  etc. 

In  general,  Bq.  (28)  cannot  be  solved  analytically,  so  that  we  seek  an  approximate  solu¬ 
tion.  First,  we  obtain  the  boaogeneous  solution  by  applying  the  equivalent  linearisation 
procedure.3  Thus,  we  have 

fo(t)  ■  cos[g(a)]  g  and  g(a)  •  g(u)t  ♦  g,  (29) 

where  q  la  the  aaplltude  vector  of  the  hath  oscillators,  and  Q  the  frequency  satrlz  and  5 
the  phase  sat r lx.  Bence,  we  can  write  down  a  first-order  formal  solution  for  s(t)  as 

g(t)  -  s^t)  ♦  njft)  ,  (30) 

where 


*l(t)  -  cos[£(a)]  q  v^dt*  g"1  sln[g(t-t,)]|QP  j(t')  (31) 

and 

*2(t)  -  j^dt*  sin[Q(t-t' )][j£(t' )  ipg  +  sfl(t)  bpq  y'(t)  ♦•••]y(f),  (32) 

and  again  sq  h  ■  a.  i.  Subsltutlng  (30)  Into  (27),  we  obtain  the  generalised  Langevln 
equation 

£(t>  -  -|pp  j(t>  -  fcpq  ^(t)  ♦  «pp  ♦  J'  fcppj'  +  •••lx(t> 

♦  J'(t)  &q  8(t)  ♦  £'  £pq  |(t)  *(t)  ♦  •••  (33) 

The  first  two  terms  on  the  BBS,  which  represent  the  harmonic  approximation,  are  equivalent 
to  the  formalism  and  Adelman  and  Doll.*  The  third  term  describes  the  enharmonic  coupling 
among  the  y-osclllators.  The  rest  are  higher-order  response  terms  which  Is  general  can  be 
neglected  If  the  set  of  y  oscillators  Is  sufficiently  big.  Thus,  Sq.  (33)  becomes 

£(t)  •  [-^p  ♦  j'  fpp  ♦  j'  £ppg'  ♦  •••lj(t)  -  |pq  s^(t).  (34) 


Bow,  ms  Include  the  motion  of  the  adspecles  to  form  a  complete  set  of  equations  of 
notion.  list  (x(t)}  be  the  coordinates  of  the  adspecles.  Be  then  have 

|<t)  -  lx  (x(t ),y(t),s)  (35) 

£(t)  •  (-|pp  ♦  X*  |PP  +  (PP  J'  ♦*••]£<*>  ♦  Iy(*(t),y(t),I)  ♦  B(t)  -£dt*  |(t-t' )j(t' )(38) 

with  fs  (x(t),y(t),7)  •  -ljrT(s(t),y(t),s)  and  |^x(t),y(t),I)-  -fyF(x(t),y(t),i). 

1  la  the  squlllhrlai  position  of  the  bath  osoillators.  Bore,  $(t)  -  knoO*1  sln(flt )hp  Is 
the  r sop os si  fusetlss  aad  B(t)  ■  -too  «os(St*|)|  Is  tbs  fluctuating  forms.  These  quantities 
are  related  by  tbs  sussed  fleetest Zod— dies 1 pit ioc  theorem* 


(37) 


<g(t)  R+(0)>  -  |pq  COS(gt)<*(0)  J+(0)>jjOP  -  kgT  g(t) 

where  we  have  used 

g(0)  ■  cos(£)a  (38) 

sad 

<j(0)  f+(0)>  -  g“al%T  (30) 

for  the  emaoaiesl  ensemble  average.  The  above  formalism  is  Ideal  for  the  study  of  the 
dynamics  of  flash  desorption  by  OV  radiation,  provided  the  excited  electronic  potential 
V(x(t ) ,y(t) ,s)  is  known.  The  initial  vibrational  temperature  of  the  y -oscillators ,  which  is 
in  general  higher  than  the  bath  temperature,  can  be  chosen  by  means  of  the  Franck-Condon 
approximation.  Bence,  one  can  study  how  these  y -oscillators  relax  and  transfer  energy  to 
the  bath  and  as  well  as  to  the  adspecies,  which  leads  to  (thermal)  desorption.  However,  ex¬ 
cited  states  associated  with  the  adspecies— solid  interaction  are  rarely  known,  which  is  the 
worst  obstacle  to  solving  the  overall  problem. 

The  above  formalism  can  also  be  used  to  study  IB  laser— induced  desorption  by  slsply  in¬ 
cluding  the  oscillating  force  (representing  the  driving  force  via  the  IB  laser)  in  both 
Bos.  (35)  and  (30).  Ve  would  like  to  point  out  that  previous  studies6**8  only  the  adbond  or 
phe  internal  motion  of  adspecies  are  IB  active.  In  another  words,  the  oscillating  force  is 
included  only  in  IQ.  (35).  Here,  we  assume  some  of  the  surface  modes  (y-oseillators)  are 
also  IB  active,  having  frequency  compatible  with  the  IB  frequency.  Actually,  this  is  a 
reasonable  assumption,  since  these  surface  atoms  which  surround  the  adspecies  behave  quite 
differently  than  the  bulk  atoms.  Thus,  the  (Ueorption  mechanism  considered  here  is  not  .lust 
a  resonance  process  as  considered  previously,7*6  but  also  an  incoherent  process  which  is  in¬ 
duced  by  the  local  heating  of  the  y -oscillators.  list  T-  be  the  steady-state  temperature  of 
the  y-oaelllators  and  On  be  the  dissociation  energy  of  the  adbond.  If  k&T.  >  Dp,  the  de¬ 
sorption  mechanism  is  as  incoherent  process.  If  kgT.  <  Do  <  kgTy  ♦  where  sy  is  the 
frequency  of  the  laser,  the  desorption  mechanism  is  then  a'  partially  coherent  process ,  as 
hhown  in  Fig.  3.  The  steady-state  temperature  of  the  y-oseillators,  which  is  in  general 
higher  than  the  bath  temperature,  is  due  to  two  complicated  completing  procesaest  multi¬ 
photon  excitation  and  multiphonon  relaxation  of  the  surface  nodes.  The  steady-state 
temperature  must  satisfy  one  of  the  two  above  conditions!  otherwise,  an  unrealistic  high- 
power  laser  is  required  to  desorb  the  adspecies.7*6  One  moat  search  for  an  optima  laser  ( 
frequency  for  the  system  of  interest  in  order  to  achieve  the  desired  steady-state 
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temperature  with  a  laser  of  reasonable  power  density.  The  steady-state  temperature  Is 
given  byt 

«,(t0) 

Hkr "° 

with 

i  fto+r 

Ty<t0)  -  jjjjj  dt  <y(t)y(t)>  . 


Here,  T Is  the  period  of  the  IB  laser. 


(40) 


(41) 


The  usual  shortcomings  arise  In  this  classical  approach.  For  example,  spontaneous 
emission  and  tunneling  effects  are  not  treated.  Fortunately,  the  deterministic  mechanisms 
of  the  above  desorption  procsss  are  the  Incoherent  heating  and  single-photon  absorption. 
Therefore,  the  neglect  of  quantum  effects  will  only  lead  to  a  small  error  in  the  calculation. 


Conclusions 

Both  quantum  mechanical  and  classical  treatments  of  the  dynamical  role  of  phonons  In  the 
excitation  and  relaxation  of  adpsecles  have  been  presented.  Using  a  linear-chain  model, 
the  former  treatment  in  Section  XI  points  out  the  Inadequacy  of  the  Markovian  approximation. 
Furthermore,  the  laser  pumping  can  be  comparable  to  multiphonon  relaxation  for  conditions 
such  as  a  shallow  potential  well  for  the  adbond  and/or  very  close  resonance  between  the 
laser  and  adbond  frequencies.  The  latter  classical  treatment  la  Section  III  Is  better  suit¬ 
ed  to  a  three-dimensional  analysis.  For  example,  here  the  response  function,  which  Is  the 
classical  analog  of  the  kernel  function  used  In  Section  IX,  is  more  easily  obtained.9 
It  is  seen  that  local  heating  of  the  surface  atom  by  an  IK  laser  is  an  Important  mechanism 
for  desorption. 
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